This report documents the method for obtaining the total organic carbon content (TOC) in the Zostera noltii seagrass. This includes information compiled from:

All work belong to the REWRITE project.

1 Method

1.1 Study sites and planning of collection points

Fieldwork was done in Cadiz, Aveiro and Bourgneuf Bay. These three sites have been covered in Davies et al. (2024).

1.1.1 Cadiz Bay

Cadiz Bay is one of the demonstrator sites for the REWRITE PROJECT. Its seagrass meadow was estimated to be 2.09 km2 during the peak season of 2024 (Davies et al. 2024).

As one of the demonstrator sites, Cadiz Bay has a rewilding/restoration site and a natural/control site. However, seagrass meadows are only viable at the natural/control site, owing to the environmental factors at the rewilding site. As such, only the natural site is of concern for this seagrass study.

1.1.2 Ria de Aveiro

Ria de Aveiro is one of the demonstrator sites for the REWRITE PROJECT. Its seagrass meadow was estimated to be 3.52 km2 during the peak season of 2024 (Davies et al. 2024).

Similar to Cadiz Bay, Ria de Aveiro is also a demonstrator site with a rewilding and natural site. Seagrasses are found in both sites. For each site, a grid sampling system was implemented.

At the rewilded site, a four by eight grid was set up, with a 50 m distance from the next station. In addition, a random sampling of 21 stations was added to the same grid area. Meanwhile, a four by seven grid was used instead, with a 50 m distance from the next station. Similar to the rewilded site, a random sampling of 18 stations was added to the same grid area.

During the fieldwork, it was discovered that some of the planned grid stations were marshes or creeks, and thus unsuitable habitats for seagrass monitoring. Thus, a total of 85 stations were sampled at Ria de Aveiro. Table 1.1 shows the breakdown by site and sampling system.

Table 1.1: Field sampling design at Ria de Aveiro. All stations at Ria de Aveiro were given a two letter code corresponding to their site and sampling system, and a station number.
Code Site Sampling system Planned no. of field stations Actual no. of field stations
RG Rewilded Grid 32 25
RR Rewilded Random 21 17
NG Natural Grid 28 27
NR Natural Random 18 16

1.1.3 Bourgneuf Bay

The seagrass meadow in Bourgneuf Bay is located close to ISOMer, the researchers’ workplace, and has been covered in multiple studies Davies et al. (2024). Its meadow was estimated to be 4.94 km2 during the peak season of 2024 (Davies et al. 2024).

A historical analysis of the meadow over 40 years was performed to identify areas in the meadows which are either consistently dense in seagrass or have a highly variable seagrass density. Site L was located in an area with a strong variability in seagrass density while site H was placed in a consistently dense part of the meadow.

For each site, 30 stations were chosen to represent a range of seagrass cover from zero to 100%, based on an in-situ visual assessment. Unfortunately, one sample for site L was unfortunately lost between the field collection and sorting, most likely during the transportation, resulting in a final sample size of 30 for site L (Table 1.2).

Table 1.2: Field sampling design at Bourgneuf Bay. All stations at Bourgneuf were given a one letter code corresponding to their site, accompanied with a station number.
Code Site Planned no. of field stations Actual no. of field stations
L Highly variable seagrass presence/cover 30 29
H Consistently dense 30 30

1.2 Sample collection

Field data collection

SPC

Hyperspectral measurement

Coring

DGPS

The collected samples were rinsed and sieved using a 40 cm diameter stainless steel sieve with a mesh size of 1 mm to remove as much substrate as possible. The samples were then transported to back to the laboratory and preserved in a fridge.

1.2.1 Sorting

1.2.2 Drying

The samples were dried in an oven at 60 °C. [ask Alex]

1.2.3 Grinding

The dried samples were placed in 50 ml falcon tubes with a ceramic ball of diameter 20 mm. The tubes were loaded onto the Mixer Mill MM 400 and grinded at 30 Hz at 30 seconds. If the initial grinding was insufficient, it was repeated at the same intensity and period until the samples were homogeneous. Homogeneous powdered samples are required for the carbon content measurement and will be elaborated in section 1.4.

1.3 Data selection

For the year 2025, a total of 161 field samples were collected across the three sites (Table 1.3). In principle, it would be best to measure the carbon content of all the available biomass samples. However, not only does each sample consist a leaf and root sample, but also both the total carbon (TC) and inorganic carbon (IC) content needs to be measured to obtain the total organic carbon (TOC) content (see section 1.4). Thus, the actual amount of laboratory work is approximately quadrupled per station with a non-zero amount of seagrass biomass. To optimise the laboratory work to a more feasible load, 30 samples each were selected from the Aveiro and Bourgneuf Bay locations, while all 17 samples in Cadiz were kept.

Table 1.3: Number of biomass samples collected from the field. Each sample consists leaf and rhizome materials, which correspond to the aboveground and belowground biomass. To ensure the representation of the range of seagrass cover, stations with zero seagrass cover are also included. Naturally, most of these stations would either have no seagrass biomass or a very low amount.
No. of field samples No. of field samples with biomass
Cadiz Bay 17 17
Ria de Aveiro 85 67
Bourgneuf Bay 59 54

1.3.1 Ria de Aveiro

The samples at Ria de Aveiro were selected based on the following factors:

  • Deviation from a linear relationship between the collected leaf and root biomass
  • Clustering in the 2D plot space between the collected leaf and root biomass
  • Deviation from a linear relationship between the seagrass cover and the leaf biomass
  • Availability of in situ hyperspectral measurements (n = 23)
  • Elimination of stations which have no seagrass biomass (n = 17)

This means that 7 samples need to be selected from the remaining 45 samples which have neither hyperspectral measurements nor have no seagrass biomass.

Figure 1.1: Data filtering for the Ria de Aveiro samples. The light blue rectangle denotes that the station has insufficient biomass for the carbon measurement, based on the requirements outlined in section 1.4. (A) establishes the general relationship between leaf and root biomass. (B) scores the non-zero stations based on the established criteria. (C) demonstrates that the filtered data’s relationship between leaf and root biomass is somewhat preserved with its similar slope factor to A (represented by the blue dotted line).

Based on the criteria, the final list of samples from Ria de Aveiro for carbon measurement are:

  • NG 10, NG 12, NG 13, NG 15, NG 16, NG 17, NG 18, NG 19, NG 23, NG 25, NG 26, NG 29, NG 6, NR 40, NR 43, NR 44, NR 46, NR 49, NR 52, RG 12, RG 13, RG 14, RG 15, RG 16, RG 22, RG 5, RG 9, RR 36, RR 38, RR 41, RR 45, RR 46.

Three stations (NG 29, NR 44, and RG 9) have too little biomass samples (Figure 1.1), and thus will be combined to ensure sufficient biomass for the carbon measurement. In addition, RG 22 and RR 46 have less than 370 mg of dry weight for either its leaf or root biomass. However, these points have a minimum of 140 mg each, which might suffice (more in section 1.4).

1.3.2 Bourgneuf Bay

The data was selected using the same criteria from section 1.3.1 for the carbon measurement. As before:

  • Availability of in situ hyperspectral measurements (n = 30)
  • Elimination of stations which have no seagrass biomass (n = 4)
  • Owing to the lack of seagrass cover data, this criterion was dropped

Normally, this would means that 0 samples need to be selected from the remaining 23 samples which have neither hyperspectral measurements nor have no seagrass biomass. Furthermore, there is a loss of 1 more sample(s), as its root biomass was unfortunately lost during the sorting and drying process. However, 7 samples have less than 370 mg for either the leaf or the root biomass or both. Of these, 1 sample(s) have less than 140 mg. This means that these samples might be pooled, so more samples are needed for the carbon measurement.

Figure 1.2: Data filtering for the Bourgneuf Bay samples. The light blue rectangle denotes that the station has insufficient biomass for the carbon measurement, based on the requirements outlined in section 1.4. (A) establishes the general relationship between leaf and root biomass. (B) scores the non-zero stations based on the established criteria. (C) shows the change in the relationship between leaf and root biomass between the original data (represented by the blue dotted line) and the filtered data.

Based on the criteria, the final list of samples from Bourgneuf Bay for carbon measurement are:

  • H_01, H_02, H_03, H_04, H_05, H_06, H_07, H_08, H_09, H_10, H_11, H_12, H_13, H_14, H_15, H_16, H_17, H_18, H_19, H_20, H_21, H_22, H_23, H_24, H_25, H_26, H_27, H_28, H_29, H_30, L_01, L_11, L_14, L_16, L_18.

Need to comment more about the samples with less than 370 mg for either AGB or BGB.

1.4 TOC measurement

The measurement of total carbon (TC) and inorganic carbon (IC) content in fly ash were performed using a SSM-5000A Shimadzu Solid Sample Combusion Unit, attached to a proprietary software for immediate data logging (Figure 1.3).

To measure the TC content in a sample, about 70g of powdered biomass was combusted at 900 °C for 8 minutes. For the IC content measurement, between 200-800 g of powdered biomass was mixed with 9ml of phosphorous acid (HPO3) and combusted at 300 °C until no carbon dioxide is detected. The increased weight of sample used was to increase the IC content measurement, as about 200 g of sample produced around 2 mg of carbon while 700 g of sample yielded around 6 mg of carbon. The measurements were scaled to the total weight of the samples using Equation (1.1).

\[\begin{equation} Carbon\,concentration\,(g\,per\,C/kg) = \frac{Measured\,Carbon\,(mg) \times Dry\,Weight\,(\%) \times 1000}{Wet\,weight\,of\,essay\,(mg)} \tag{1.1} \end{equation}\]

The difference between the TC and IC content is the total organic carbon content (TOC).

User interface for the Shimadzu proprietary software. This shows the end result of an inorganic carbon (IC) content combustion for the leaf biomass sample number 5 from Cadiz. Using a sample load of 763.2 mg with 9 ml of HPO~3~, an initial IC content measurement of 6.267 mg was recorded. The actual IC content of the sample would still need to be computed using Equation \@ref(eq:carbon).

Figure 1.3: User interface for the Shimadzu proprietary software. This shows the end result of an inorganic carbon (IC) content combustion for the leaf biomass sample number 5 from Cadiz. Using a sample load of 763.2 mg with 9 ml of HPO3, an initial IC content measurement of 6.267 mg was recorded. The actual IC content of the sample would still need to be computed using Equation (1.1).

Based on our visual observation, it seemed important for the samples to be sufficiently homogeneous for two reasons. Firstly, more samples can be added to the sample boat if it were homogenised into powder rather than having bits of larger biomatter creating space within. Based on our limited laboratory experience, we found that only about 200-300 mg of incompletely homogenised sample can be placed into the SSM 5000A sample boat, while up to about 800 mg of homogenised sample can be added. Secondly, it was observed for the TC contect combustion that the insufficiently homogenised samples have more leftovers than the homogenised samples (Figure 1.4). This suggests the possibility of either a slower or incomplete combustion, which may affect the measurements.

Samples post-combustion for the total carbon (TC) content measurement. The sample on the left was not sufficiently homogeneous prior to combustion while the sample on the right was finely ground. The relatively larger pieces of biological matter on the left is still perceivable even after combustion at 900 °C. Furthermore, by visual analysis, its remaining content seems to be more than the homogeneous sample on the right. Do note that this is purely a visual observation based on the sample homogeneity pre-combustion and their corresponding post-combustion state.

Figure 1.4: Samples post-combustion for the total carbon (TC) content measurement. The sample on the left was not sufficiently homogeneous prior to combustion while the sample on the right was finely ground. The relatively larger pieces of biological matter on the left is still perceivable even after combustion at 900 °C. Furthermore, by visual analysis, its remaining content seems to be more than the homogeneous sample on the right. Do note that this is purely a visual observation based on the sample homogeneity pre-combustion and their corresponding post-combustion state.

References

Bargain, A, M Robin, V Méléder, P Rosa, E Le Menn, N Harin, and L Barillé. 2013. “Seasonal Spectral Variation of Zostera Noltii and Its Influence on Pigment-Based Vegetation Indices.” Journal of Experimental Marine Biology and Ecology 446: 86–94. https://doi.org/10.1016/j.jembe.2013.04.012.
Davies, Bede Ffinian Rowe, Simon Oiry, Philippe Rosa, Maria Laura Zoffoli, Ana I. Sousa, Oliver R. Thomas, Dan A. Smale, et al. 2024. “Intertidal Seagrass Extent from Sentinel-2 Time-Series Show Distinct Trajectories in Western Europe.” Remote Sensing of Environment 312: 114340. https://doi.org/10.1016/j.rse.2024.114340.
Zoffoli, Maria Laura, Pierre Gernez, Philippe Rosa, Anthony Le Bris, Vittorio E Brando, Anne-Laure Barillé, Nicolas Harin, et al. 2020. “Sentinel-2 Remote Sensing of Zostera Noltei-Dominated Intertidal Seagrass Meadows.” Remote Sensing of Environment 251: 112020. https://doi.org/10.1016/j.rse.2020.112020.